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It is shown that bulk ferromagnets support propagating non-linear modes that are analogous to 
the vortex rings, or "smoke rings", of fluid dynamics. These are circular loops of magnetic vorticity 
which travel at constant velocity parallel to their axis of symmetry. The topological structure 
of the continuum theory has important consequences for the properties of these magnetic vortex 
rings. One finds that there exists a sequence of magnetic vortex rings that are distinguished by 
a topological invariant (the Hopf invariant). We present analytical and numerical results for the 
energies, velocities and structures of propagating magnetic vortex rings in ferromagnetic materials. 
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In many situations ferromagnetic materials may be 
viewed as continuous media, in which the state of the 
system is represented by a vector field indicating the 
local orientation of the magnetisation. The dynamics 
of the ferromagnet then follow from the time evolution 
of this vector field, which obeys a non-linear differential 
equation known as the Landau-Lifshitz equation. Repre- 
senting the local orientation of the magnetisation by the 
vector field n(r, t) (n is a three-component unit vector), 
the Landau-Lifshitz equation [hj takes the form 
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in the absence of dissipation, p is the density of mag- 
netic moments, each of angular momentum J, and E is 
the energy, which is a functional of n(r, t) and its spatial 
derivatives (summation convention is assumed through- 
out this paper). 

The Landau-Lifshitz equation has been the subject of 
numerous studies , owing to its physical importance 
as a general description of ferromagnetic materials, and 
to the rich mathematical properties that result from its 
combination of non-linearity and non-trivial topology. Of 
particular interest are the solitons and solitary waves |t]] 
that it has been found to support. In one spatial di- 
mension, the Landau-Lifshitz equation is integrable for 
certain energy functionals and the complete set of soli- 
tons is known In higher dimensions, the equation 
is believed to be non-integrable for even simple energy 
functionals || , and the understanding of non- linear exci- 
tations g is incomplete. 

Here, we construct a novel class of solitary waves of 
three-dimensional Landau-Lifshitz ferromagnets. Our 
approach relics on the conservation of the linear momen- 
tum 



(2) 



P a = 2TrpJe a /3ry / d r rpVL 1 {r, t) 



the definition of the hydrodynamic impulse of an incom- 
pressible fluid jyj, if fluid vorticity is identified with 
magnetic vorticity. Our work emphasises the connection 
to fluid dynamics by showing that there exist solitary 
waves in ferromagnets that are analogous to the vortex 
rings of fluid dynamics jlljjl^]: these are circular loops 
of (magnetic) vorticity that propagate at a constant ve- 
locity parallel to their axis of symmetry. (The magnetic 
analogues of vortex/anti-vortex pairs have recently been 
determined using a similar approach fl3| .) 



where il a = -^eap^eijkniX/ pn 3 ;V 7 n/c is the "magnetic 
vorticity" . This definition of momentum [fill resembles 



Furthermore, just as there exist generalisations |L4 15 
of the vortex rings in fluid dynamics to vortex ring struc- 
tures in which the lines of vorticity are linked (as mea- 
sured by a non-zero helicity |l6f| ) , there exist similar gen- 
eralisations of the magnetic vortex rings in ferromagnets 
to topologically non-trivial structures involving the link- 
ing of vortex lines |lj],[l^] (in this context the measure of 
linking is known as the Hopf invariant jl9| ) . As a conse- 
quence, ferromagnets support a sequence of topologically- 
distinct magnetic vortex rings. One can understand how 
this sequence arises by noting that a magnetic vortex 
line carries an internal orientation that can be twisted 
through an integer multiple of 2n as the vortex line traces 
out a closed loop |^| . By varying the number of rotations 
of this internal angle, one obtains a sequence of magnetic 
vortex ring configurations that are topologically distinct 
(they cannot be interconverted by non-singular deforma- 
tions), as they relate to different values of the Hopf in- 
variant p9| , (See Ref. |^l| for a similar construction for 
fluids, where topology allows the inserted twists to be 
non-integer multiples of 2tt; magnetic vortex rings are 
more akin to the coreless vortex rings of superfluid 3 He- 
A, Ref. ^2|.) The Hopf invariant, H, is the integer invari- 
ant that characterises the mappings S 3 — > S 2 (the vector 
field n{r) describes such a mapping when fixed bound- 
ary conditions are imposed, e.g. n(\r\ — > oo) = +z); it 
can be interpreted in terms of the linking number of two 
vortex lines on which n takes different values |ll| . This 
topological invariant has been of interest in recent stud- 
ies of non- linear field theories, where it has been used 
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to stabilise static solitons 23-2511. Here we show that it 



classifies a sequence of dynamical solitary waves of ferro- 
magnets - the magnetic vortex rings. 

We shall construct magnetic vortex ring solitary waves 
for ferromagnetic materials described by the energy func- 
tional 
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which represents isotropic exchange interactions and uni- 
axial anisotropy (we consider only A > 0, and choose the 
groundstate to be the uniform state with n = +z). It 
is straightforward to verify that, for this functional, the 
momentum (^) is conserved by the dynamics as is 
the number of spin-reversals 



N = pJ/h / d 3 r (1 -n z ). 
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(Within a full quantum description, the number of spin- 
reversals would be an integer; within the semiclassical de- 
scription afforded by the Landau-Lifshitz equation, which 
is accurate for N 3> 1, N a continuous variable.) Our 
approach is to find configurations, n*(r), that extrem- 
ise the energy (^) at given values of the momentum, P, 
and number of spin-reversals, N, within each topologi- 
cal subspace, TL. This procedure |^6|,|l3| defines an ex- 
tremal energy E^ i (P 1 N). The variational equations can 
be used to show that there exist time-dependent solutions 
of Eqn.(Q) of the form 



n x (r,t) + in y (r,t) = [n* x (r - vt) 
n z (r,t) = n* z (r - vt) 
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These solutions describe travelling waves which move in 
space at constant velocity, v a , while the magnetisation 
precesses around the z-axis at angular frequency uo. We 
shall find configurations n*(r), resembling magnetic vor- 
tex rings, which have spatially-localised energy density 
and therefore describe propagating solitary waves . 

It has been suggested previously that solitary waves 
resembling magnetic vortex rings might exist in ferro- 
magnets, being stabilised by a non-zero Hopf invariant 
TL combined with the conservation of either the number 
of spin-reversals JlTj N, or the momentum |18j P. In 
the present work we make use of the conservation of both 
N and P, allowing solitary waves with more general mo- 
tions to be constructed. This proves to be essential for 
the magnetic vortex rings that we discuss here. Further- 
more, since we do not invoke topological stability, we can 
obtain solitary waves for all values of the Hopf invariant, 
including TL = 0. 

We now turn to the determination of the properties 
of magnetic vortex ring solitary waves within the dif- 
fering topological subspaces TL. First, we make some 



general remarks. Due to the invariance of the en- 
ergy (H) under spatial rotations, the extremal energy 
E?h(P, N) is independent of the direction of momentum; 
for later convenience, we choose P = Py. By consid- 
ering the dependences of Eqns.(|2|||||) under scale trans- 
formations, one finds that the extremal energy has the 
form E^{P,N) = N x l 3 &H{p,a), where p = P/N 2 / 3 and 
a = AN 2 / 3 (here, and subsequently, we choose units for 
which p s = pj = h = 1). 

As a first step, consider linearising the equations of 
motion about the groundstate n(r) = z. The configu- 
rations that extremise the energy at fixed N and P are 
easily found: they are spatially-extended, and describe 
N non- interacting spin- waves each of momentum P/N. 
The results of this analysis determine the spin-wave dis- 
persion Sj^Lq (p, a) = p 2 + a (the Hopf invariant is zero 
for all small-amplitude disturbances). 

Insight into the results of the full, non-linear theory 
may be achieved by considering large-radius magnetic 
vortex ring configurations. Consider a magnetic vortex 
carrying a total "flux" Q, that is closed to form a cir- 
cular loop with radius much larger than the size of the 
vortex core. (The vortex core size, like the loop radius, 
varies with N and P. The condition that the core is small 
compared to the radius is P/QN 2/3 = p/Q > 1.) The 
flux Q is defined to be the integral of the vorticity across 
a surface pierced by the vortex; since the magnetisation 
tends to a constant, n = z, far from the vortex core, Q is 
an integer (this is the topological invariant that classifies 
the mappings S 2 — ► S 2 ). Neglecting the size of the vor- 
tex core in comparison to the radius of the loop, one can 
determine the minimum energy of these magnetic vortex 
rings using results from previous studies of topological 
solitons in two-dimensional ferromagnets H . One finds 
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S^ m (p,a) = 2V2^(p/|Q|) 1 / 4 , 



a 2 <£Lp/Q, 
a 2 > p/Q, 
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which are valid for p/Q 3> 1, when the vortex core is 
small compared to the radius of the vortex ring (Q ^ 
has been assumed). The relative sizes of a 2 and p/Q 
determine whether the vortex core is small [Eqn. m)] 
or large [Eqn. (Q)] compared to a lengthscale set by the 
anisotropy (the width of a domain wall). 

In order to determine the properties of the magnetic 
vortex rings when the size of the vortex core is compa- 
rable to the radius of the ring (p/Q ~ 1) we have em- 
ployed numerical analysis (we have not been able to solve 
the variational equations analytically). We simplify the 
problem by minimising the energy within a class of ax- 
isymmetric configurations j24| that is consistent with the 
general variational equations: 

Mr) + in y (r)] = [n x {p, y) + in y {p, y)} e lM *, y > 0, 
Mr) + m„(r)] = [n x (p, -y) - in y (p, -y)] e" lM0 , y < 0, 
n z (r) = n z (p, \y\), 
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where (p, y, 4>) are cylindrical polar co-ordinates about 
the y-axis. We thereby reduce configuration space to the 
functions n(p, y) on the quarter plane (p > 0, y > 0) and 
the integer parameter M. The Hopf invariant of these 
configurations is p3 Ti = MQ, where Q is the total flux 
through the half-plane (p > 0). For M ^ 0, finite en- 
ergy configurations have n — > z at p — ► as well as at 
p 2 + y 2 — > oo, so the total flux Q, and hence the Hopf 
invariant, take only integer values. 

We have studied the isotropic ferromagnet, a — 0, by a 
discretisation of the region (0 < p < L,0 < y < L/2) on 
square lattices of up to 600 x 300 sites. Fixed boundary 
conditions, n = z, were imposed on p = L and y = L/2, 
and n y = was imposed on y — for consistency with the 
above ansatz (the requirement that n = z on p = when 
M =/= emerges naturally from the energy minimisation) . 
Energy minimisation was achieved by a conjugate gradi- 
ent method with constraints on N and P imposed by an 
augmented Lagrangian technique \27\. 
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FIG. 1. Energies of the magnetic vortex ring solitary waves 
with Ti — and TL = 1 for a — (shown as the ratio to 
the asymptotic energy ([J) with Q = 1). To indicate the ex- 
tent of finite-size effects, two system sizes are shown: the 
largest we have studied L = 11.2N 1 / 3 /p 1 ^ (solid lines), and 
L — 8.2./V 1/ ' 3 /p 1/ ' 4 (crosses). The spin-wave dispersion is 
shown as a dashed line. Inset: the same quantities as func- 
tions of l/p 3/2 . 

The energies of the magnetic vortex ring solitary waves 
with H = (M = 0) and H = 1 (M = Q = 1) are shown 
in Fig. ^ as functions of the scaled momentum p. At 
large momenta, both branches approach the asymptotic 
expression (^|) for a Q = 1 vortex. Assuming that the 
leading corrections shown in the inset are linear in 1/p 3 / 2 , 
one finds from Eqns. (|^,^) that a vortex loop of large ra- 
dius, p/Q 3> 1, moves in such a way that its precession 
frequency to is comparable to the frequency with which 
it translates a distance of order its own radius. At small 
momenta, the solitary waves become higher in energy 
than non-interacting spin-waves for p < 4.66 (Ti = 0) 



and p < 5.1 (Ti = 1). Both branches are found to persist 
below these values for a small range of p as local energy 
minima. 

Typical configurations of the solitary waves are illus- 
trated in Figures || and ^, which show the local magneti- 
sation within the xy-plane and three-dimensional repre- 
sentations of the curves on which n = — z and n = +x. 
These curves allow a direct visualisation of the topolog- 
ical structure of the configurations, as measured by the 
Hopf invariant |Tjfl . In Fig. ^|, the curves are simply cir- 
cles centred on the ?/-axis: they are unlinked, illustrating 
that the Hopf invariant is zero, Ti = 0. In Fig. ||, the two 
curves are linked once, illustrating a non-trivial topolog- 
ical configuration with unit Hopf invariant, Ti = 1. 




FIG. 2. The configuration of the magnetic vortex ring with 
TL = 0, p = 5.0, and a = 0. The vectors represent the projec- 
tion of n onto the xy-plane; their colours indicate the value of 
n z (refer to the colour- legend). The tubular curves illustrate 
the lines along which n = — z (red) and n — +x (green). (For 
clarity, only the central section of the full lattice is illustrated, 
with a fraction of the lattice points reproduced.) 




The same as Fig.^, for Ti = 1, p = 5.3, and a 



For larger values of the Hopf invariant (Ti > 2) we find 
that corrections arising from the finite lattice size in our 
calculations become more significant. (Finite-size effects 
are already apparent in Fig. [I] for Ti = 1.) These effects 
prevent a convincing demonstration of the existence of 
(non-singular) magnetic vortex rings with Ti > 2 that 
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describe stable energy minima. The construction of sta- 
ble magnetic vortex rings with Ti. > 2 may require the 
use of non-axisymmetric configurations. The results of 
the present study demonstrate that magnetic vortex rings 
with H = and 7i = 1 do describe stable energy minima 
within the class of axisymmetric configurations assumed. 
Our results provide the structures, energies, and there- 
fore the velocities and precession frequencies [see Eqns. 

] of these propagating non-linear modes. Further nu- 
merical studies p8| indicate that similar magnetic vortex 
ring solitary waves exist for non-zero anisotropy, a =/= 0. 
(In fact, the magnetic vortex rings are apparently more 
favourable: e.g. the branch of solitary waves with H = 
persists down to p = for a > 23.5, consistent with the 
existence of purely precessional solitary waves in a uniax- 
ial ferromagnet jq|.) Any additional sources of magnetic 
anisotropy in experimental systems, or the inclusion of 
magnetic dipole interactions |2^] , will lead to corrections 
that are small when the magnetic vortex rings are suf- 
ficiently small (compared to a characteristic lengthscale 
set by the strength of these additional couplings) . 

One may wonder why magnetic vortex rings have not 
as yet been observed experimentally, whilst vortex rings 
in fluids are a matter of everyday experience. The an- 
swer lies in the difficulty of creating non-linear excitations 
in solid-state materials. For example, magnetic vortex 
rings involve a large number of spin-reversals, so they 
are not accessed in standard inelastic neutron scatter- 
ing experiments which probe single spin- flip excitations. 
The creation of magnetic vortex rings in ferromagnetic 
materials will require the use of other experimental tech- 
niques. One way in which magnetic vortex rings could be 
excited experimentally, the details of which we are cur- 
rently investigating, is to exploit an instability that we 
have discovered to the creation of magnetic vortex rings 
in itinerant ferromagnets of mesoscopic size under con- 
ditions of high current density. This instability signals a 
transition to a form of magnetic turbulence in mesoscopic 
ferromagnets, driven by the exchange coupling between 
the magnetisation and the spins of the itinerant electrons, 
and may be relevant to the unexplained dissipative phe- 
nomena observed in Ref. (3^] . 

The author is grateful to Mike Gunn and Richard Battye 
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